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Disturbances Due to the Motion of a Cylinder in a 

Two-Layer Liquid System 1 
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The disturbance created at the interface of a two-layer liquid system by the horizontal 
motion of a cylinder in the upper layer is studied for various sizes and shapes of cylinders, 
depths of the liquids, cylinder velocities, and density ratios. The disturbances fall into 
three categories. First, when the layers are of equal thickness, in most cases a train of 
progressive oscillatory waves is produced at the interface. Second, when the depth of the 
denser layer is much less than the depth of the fresh-water layer, the profile of the interface 
usually resembles that corresponding to positive internal solitary waves. Third, when the 
depth of the denser layer is much greater than the depth of the fresh-water layer, in most 
cases an internal hydraulic jump is produced. The characteristics of the disturbances 
in each category are related to the size of the cylinder, the depths of the liquids, the cylinder 
velocity, the density ratio, and the total distance of travel of the cylinder. A theoretical 
analysis is given for disturbances of the first category. 



1. Introduction 

Almost a century ago Stokes gave the theory of 
internal waves of infinitesimal height in a two-layer 
system, the lighter liquid resting on the heavier one 
[l], 2 and about a quarter century later the problems 
of multilayer systems were considered by Webb 
[2], Grenhill [3], and Burnside [4]. Following these, 
the wave motions of hetrogeneous liquids with 
continuously varying density gradients in the vertical 
were taken up by Love [5], "Rayleigh [6], and Lamb 
[7], These studies relate to an endless succession of 
waves, standing or progressive, without reference to 
the generation of the disturbance such as would result 
from the motion of objects through a stratified 
medium or from the uniform flow of fluids over a 
ridge fixed on the bottom. 

The first studies of the latter category go back 
to Ekman, who considered disturbances of the inter- 
face of a two-layer system initiated when a body is 
towed over the surface [8]. The occasion of the 
study was the classical account by Nansen of dead 
water experienced by the Fram during the North 
Polor Expedition near the coast of Siberia. Ekman 
approached the problen with much encouragement 
from Bjerknes both from an experimental and 
theoretical point of view. In addition to consti- 
tuting the first serious attempt at model procedures 
in stratified fluids, the investigation is significant 
in two respects. First, working with a model of 
the Fram, Ekman showed that the internal waves 
brought into play an added resistance having a 
magnitude dependent on the ratio of draft to fresh 
water depth, the ratio of layer depths, and the pro- 
portionate density differences. The maximum 
values of the added resistances were obtained when 
the ratio of the velocity of the model to the velocity 
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of long internal waves of the bilayer system attained 
the value 0.73, and the added resistance gradually 
decreased for larger velocities of the model. Second, 
the question of the added resistance was considered 
using a method of approach originally adopted by 
Kelvin in a similar problem applied to a homogeneous 
fluid. In the evaluation the complete representation 
of the rufflement of the interface was avoided, and 
attention was focused only on the parts contributing 
to the resistance. In this theory the incidence of 
the maximum is correctly predicted. 

The converse problem, that is the problem of a 
two-layer system moving over a ridge, was considered 
by Zeilon [9]. The investigation had its inspiration 
in the views of Pettersson to the effect that the 
internal boundary waves resulting from the interac- 
tion of the Atlantic tides on the irregular bottom con- 
figurations had a dominating influence upon the 
circulation of the Baltic. The problem investigated, 
thus, was the perturbed state of the interface owing 
to a ridge at the bottom of a two-layer system in 
simple harmonic horizontal motion. Experiments 
were carried out which were intended to verify in 
a qualitative sense the results of the lengthy mathe- 
matical analysis. 

Among the recent studies the contributions of 
Robert Long have much to offer. A novel aspect is 
the first integration of the hydrodynamic equation 
for the two-dimensional case of nonviscous flow con- 
fined to the space between two rigid horizontal 
boundaries, the initial velocities (always horizontal) 
and the densities of the undisturbed regions being 
prescribed. Rigorous solutions were obtained for 
the cases where the liquid flows initially with uni- 
form velocity and with linearly varying densities 
over a flat rectangular ridge [10]. The analysis 
was repeated for a smooth object, and experiments 
were made which showed a wave motion giving 
qualitative confirmation of the theoretical results 
[11]. The interpretations of the theory and ex- 
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periments have numerous aspects, but it should suf - 
fice to state here that in general the criterion of the 
densimetric Froude number (to be defined later) 
equaling 1/t divides the types of internal flow. Of 
particular interest is the examination, experimental 
and theoretical, of the flow of two superimposed 
layers over a ridge. Internal waves were generated 
when the ridge was of small projection, whereas 
with the much enlarged and elevated ridges the most 
significant event was the formation of leeward hy- 
draulic jumps [12]. 

The question of critical flows and the manifes- 
tation of hydraulic jumps in a bilayer or a multi- 
layer system is an important one. The significance 
of the occurrence to meteorology was first empha- 
sized by Tepper [13]. The matter has been ap- 
proached from various points of view by Rossby [14], 
Craya [15], Benton [16], Yih [17], and others. 

The matter of finite amplitude waves for a two- 
layer system has been considered by Keulegan and 
Long [18, 19]. Theory predicts the existence of 
positive internal solitary waves for the condition 
that the upper layer has the greater depth. The ex- 
periments by Keulegan tend to confirm the results 
of theory relating to shape and dependence of ve- 
locity of propagation on wave height [20]. 

In common with the disturbances produced in a 
nonhomogeneous system by liquid moving over 
ridges on the bottom or of objects being towed on the 
surface, other related phenomena would be associated 
with bodies moving through the interior of the 
medium. The simplest situation arises with the 
two-dimensional flow of a two-layer system, and this 
is the subject of the present study. We are here 
concerned with disturbances of the interface pro- 
duced by a cylinder moving with constant hori 
zontal velocity in the upper layer as illustrated in 
figure 1. The liquids are initially at rest, and the 
density difference is small. Experiments are made 
with various depths of the layers, sizes and shapes of 
cylinders, and cylinder velocities. Further experi- 
ments with continuous density gradients will be 
described in a later report. 

2. Apparatus and Procedures 

The experimental channel is 26.25 m long, 23.3 cm 
wide, 57 cm deep, and is made up of 21 glass-paneled 



sections, each 125 cm long. The salt water was 
mixed and dyed with acid chrome blue and then ad- 
mitted to the channel to a depth of 2 or 3 cm above 
the desired level of the interface. The salt water 
was circulated in the channel by a circulating pump 
to assure uniform density of the salt layer. After 
the mixing was complete and the circulation had 
subsided, fresh water for the upper layer was intro- 
duced at the center of the channel through the de- 
vice shown in figure 2. The device was lowered 
until the lips touched the water surface. The fresh 
water was introduced very slowly at first so as to 
produce as little mixing as possible. A siphon was 
placed at each end of the channel to remove the 
mixed liquid at the interface. The process was con- 
tinued until the interface had been reduced to the 
desired depth, h, and the total depth had increased 
to the desired value, h-\-h' (see fig. 5). 

After the channel was filled and all internal mo- 
tion of the liquid subsided, the density distribution 
was determined by traversing vertically with a salin- 
ity electrode. Various shunts were used in the cir- 
cuit for greater accuracy in determining salinities 
within various ranges. The electrode was cali- 
brated with each shunt so the milliampere current 
readings could be used to determine the density of 
salt solutions. An example of the vertical distri- 
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Figure 2. Longitudinal section of device for introducing fresh 
water on top of saline water. 




Figure 1. Disturbance at the interface between two liquids produced by a cylinder moving with constant velocity in the upper layer. 

The symbol T indicates the water surface. 
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bution of density is given in figure 3. The thick- 
ness of the layer of mixed liquid was 1 to 1.5 cm in 
all cases which was small compared to the depths of 
the fresh and salt layers. 

(Minders of three different diameters, 2.54, 5.08, 
and 10.16 cm, and three different shapes, circular, 
lenticular, and tear drop, were used. The lenticular 
and tear drop shapes are shown in figure 4. In 
all cases the center of the cylinder was placed 30.0 
cm from the bottom. 

The cylinders were towed by an electric motor 
through a variable speed transmission, a system of 
pulleys, and two 0.81 -mm phosphor bronze towing 
wires. Four 1.22-mm piano wires under tension 
were used as guide wires to prevent vertical oscilla- 
tions of the cylinder. Through the variable speed 
transmission any velocity from to over 1 m/sec 
could be selected. The cylinder velocity was de- 



0.50* 


=0— O 


- — 




o 


OQ) 
O 


0.25 



































APb 



Figure 3. Example of vertical distribution of density. 

y is the vertical coordinate with origin at the interface. Ap b is the difference 
between the densities at the bottom and at the free surface. 




Figure 4. Cross sections of lenticular and tear drop cylinders. 



termined by timing the (ravel over a distance of 
5 m, the first reference mark being 2 m from the 
stalling position. 

The profile of the interface during the run was 
recorded with three fixed cameras which con Id he 
triggered simultaneously or individually. For the 
majority of the runs, the pictures were la ken after 
the cylinder had traveled about 8 m. 

3. Scope of the Experiments 

In an experimental study relating to the effects 
of a cylinder when moving in one of the layers of a 
two-layer system, the items to consider first are the 
state of rufflement of the interface, the deformation 
of the free surface, and the particle velocities in 
the parts far removed from the cylinder. The 
examination and the measurement of the last two 
quantities may be omitted in a first study, as these 
may be inferred from the interface deformation on 
the basis that the flows are irrotational and the 
liquids incompressible. A second item would be 
the additional resistance due to the effect of the 
internal waves. The observation of this quantity 
may be omitted, likewise, as it is amenable to a 
close evaluation from a description of the interface, 
using for example, the method of Ekman [8]. 

Accordingly, in the present study only the interface 
deformation is considered (see fig. 5). The two 
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Figure 5. Notation diagram for inlerfacial disturbance pro- 
duced by a cylinder moving with constant horizontal velocity 
in the upper layer. 

fundamental quantities to be established were the 
wave height and wave length in relation to the 
cylinder speed, cylinder position, the thickness of 
layers, and the proportionate density difference. It 
was very surprising to find that the tests in their 
entirety failed to show a unique type of disturbance. 
In effect, three modes of deformation could be 
distinguished. If the layers were nearly of the 
same thickness, then the interface in most cases 
yielded a picture of numerous progressive oscillatory 
waves. If the lower layer was of a smaller thickness 
than that of the upper, the cylinder was usually 
followed by an internal solitary wave. On the 
other hand, if the upper layer was much thinner 
than the lower, the interface ordinarily deformed 
into a hydraulic jump. 

4. Internal Oscillatory Progressive Waves 

For those runs of the study for which the salt 
water and fresh water layers were of equal depth, 
h' = h, progressive waves were produced at the 
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interface as illustrated in figure 1. The wavelength, 
X, and the semiwave height, a, of the internal waves 
depend on the following variables: 

U, velocity of cylinder; 

r, radius of cylinder; 

h, depth of saline layer; 

h' ', depth of fresh layer; 

P, density of saline layer; 

p ', density of fresh layer; 

g, acceleration of gravity; 

s, vertical distance of center of cylinder above the 

undisburbed level of the interface ; 
/J, shape parameter of cylinder; 
L, distance traveled by cylinder at time of 

observation. 

Viscous effects are not considered . For the densities 
used, p'=l and 1.02 <p< 1.08, we may write 

p—p' _Ap_Ap 

P + P P P 

The values of the variables for the various runs 
of this class, h=h' , are given in table 1. Now 

X=\ (u, r, h, ^ g, s, 0, i) 



or 



i j 






In figures 6, 7, and 8 plots are given of \/h against 
U 2 /(Ap/p)gh, a form of densimetric or internal Froude 
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Figure 6. Variation of the wavelength of internal progressive 
waves, t=1.27 cm. 



number. The wavelengths for the experiments 
were taken as the mean value for three or four waves 
immediately following the cylinder. The line repre- 
senting the theoretical wavelength is taken from 
Lamb, Art. 231 eq (17) [21], which represents wave 
motion of two superposed fluids with the upper 
surface free. For the case 

P 
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Figure 7. Variation of the wavelength of internal progressive 
waves, y = 2.54 cm. 
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Figure 8. Variation of the wavelength of internal progressive 
waves, i = 5.08 cm. 
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Table 1. Variables for h'=h 



Run 


Cylin- 
der a 


h' = h 


S 


Ap 
P 


U 


a 


X 






cm 


cm 




cm/sec 


cm 


cm 


1 


2.54A 


25 


5 


0. 02 


HA 


3.9 


40.0 


2 


2.54 


25 


5 


.02 


12. 1 


5.3 


101.0 


4 


2.54 


25 


5 


.02 


15.0 


2.5 




5 


2.54 


25 


5 


.02 


4.0 


1.2 


~~ - 10.~5 


10 


2.54 


25 


5 


.02 


9. 1 


3.6 


48.5 


11 


2.54 


25 


5 


.02 


10.0 


4.0 


70. 


12 


2.54 


25 


5 


.02 


13.9 


2.3 


115.0 


15 


2.54 


25 


5 


.04 


12.0 


3.5 


43.5 


16 


2.54 


25 


5 


.04 


13.8 


4.3 


65. 5 


17 


2.54 


25 


5 


.04 


15.7 


4.7 


87.5 


18 


2.54 


25 


5 


.04 


18.0 


5.2 


107.0 


19 
21 


2.54 
2.54 


25 


5 


.04 


21.3 


2.2 




25 


5 


.08 


16.9 


4.3 


44.3 


22 


2.54 


25 


5 


.08 


18.9 


3.5 


59.5 


23 


2.54 


25 


5 


.08 


23.9 


5.1 


102. 5 


24 


2.54 


25 


5 


.08 


22.2 


5.0 


90.5 


25 
71 


2.54 
2.54 


25 
26 


5 


.08 


27.6 


2.2 




4 


.04 


13.7 


4.4 


"~~52.~7 


72 


2.54 


26 


4 


.04 


19.2 


3.3 


111.0 


75 


2.54 


26 


4 


.04 


12.8 


4.2 


46.3 


76 


2.54 


26 


4 


.04 


17.0 


4.5 


97.5 


77 


2.54 


26 


4 


.04 


22.7 


1.3 





78 


2.54 


24 


6 


.04 


12.8 


3.9 


50. 6 


79 


2.54 


24 


6 


.04 


10.2 


3.3 


26. 2 


80 


2.54 


24 


6 


.04 


11.3 


3.5 


36. 3 


81 


2. 54A 


24 


6 


.04 


14.1 


4.1 


66. 1 


82 


2.54 


24 


6 


.04 


15.7 


5.0 


87.5 


83 


2.54 


24 


6 


.04 


19.5 


2.5 


107.0 


85 


2.54 


23 


7 


.04 


10.4 


2.8 


29.3 


86 


2.54 


23 


7 


.04 


11.3 


3.3 


36.0 


87 


2.54 


23 


7 


.04 


12.8 


3.8 


53. 


88 


2.54 


23 


7 


.04 


13.8 


3.8 


61.3 


89 


2.54 


23 


7 


.04 


16.2 


4.8 


89.0 


90 


2.54 


23 


7 


.04 


19. 1 


2.8 


107.0 


91 


2.54 


23 


7 


.04 


21.6 


1.4 




92 


2.54 


23 


7 


.04 


9.5 


1.8 


26. 4 


93 


2.54 


23 


7 


.04 


5.7 


0.4 


li. 1 


94 


2.54 


23 


7 


.04 


6. 7 


1.2 


14.8 


95 


2. 54 


23 


7 


.04 


7.8 


1.7 


16.7 


96 


2.54 


20 


10 


.04 


6.4 


0.6 


14.4 


97 


2.54 


20 


10 


.04 


7.7 


1.1 


19.8 


98 


2.54 


20 


10 


.04 


9.8 


2.1 


28.8 


99 


2.54 


20 


10 


.04 


11.3 


2.7 


39.3 


100 


2.54 


20 


10 


.04 


13.0 


3.4 


52.0 


101 


2.54 


20 


10 


.04 


13.8 


3.8 


61.0 


102 


2.54 


20 


10 


.04 


16.7 


4.0 


98.7 


103 


2.54 


20 


10 


.04 


19.2 


2.2 


105.0 


104 


2.54 


18 


12 


.04 


6.7 


0.5 


L6.6 


105 


2.54 


18 


12 


.04 


7.8 


1.0 


22. 2 


106 


2.54 


18 


12 


.04 


10.0 


1.8 


33.8 


107 


2.54 


18 


12 


.04 


11.4 


2.7 


44.6 


108 


2. 54A 


18 


12 


.04 


12.4 


3.2 


53.5 


109 


2.54 


18 


12 


.04 


14.1 


3.4 


68.3 


110 


2.54 


18 


12 


.04 


15.6 


3.8 


80.5 


112 


2. 54 


16.25 


13.75 


.04 


7.0 


0.7 


19.8 


113 


2.54 


16.25 


13.75 


.04 


9.5 


2.8 


32.2 


114 


2.54 


16. 25 


13.75 


.04 


10.3 


3.2 


37.0 


115 


2.54 


16. 25 


13.75 


.04 


11.6 


5.2 


49.4 


116 


2.54 


16.25 


13. 75 


.04 


13.0 


6.4 


63.7 


117 


2.54 


16.25 


13.75 


.04 


14.2 


6.3 


71.3 


118 


2.54 


16. 25 


13.75 


.04 


16.6 


4.1 


84.0 


142 


2. 54B 


25 


5 


.04 


11.1 


3.4 


35.5 


143 


2.54 


25 


5 


.04 


16.6 


5.7 


89.3 


144 


2.54 


25 


5 


.04 


20.0 


3.2 


120.0 


145 


2.54 


20 


10 


.04 


7.7 


1.0 


19.2 


146 


2.54 


20 


10 


.04 


13.7 


4.5 


66.7 


147 


2.54 


20 


10 


.04 


19.7 


1.7 


102.0 


148 


2.54 


16.25 


13.75 


.04 


7.8 


1.0 


24.0 


149 


2.54 


16.25 


13.75 


.04 


12.6 


4.6 


61.5 


150 


2.54 


16.25 


13.75 


.04 


15.5 


4.7 


82.2 


153 


2.54C 


25 


5 


.04 


11.1 


3.5 


34.8 


154 


2.54 


25 


5 


.04 


16.4 


4.2 


96.0 


155 


2.54 


25 


5 


.04 


20.2 


2.7 


107.0 


156 


2.54 


20 


10 


.04 


7.9 


0.8 


21.1 


157 


2.54 


20 


10 


.04 


13.9 


3.8 


73.7 



Run 


Cylin- 
der » 


h' = h 


s 


Ao 
P 


U 


a 


X 






cm 


cm 




cm/sec 


cm 


cm 


158 


2.54 


20 


10 


AH 


19.3 


2.1 


97.5 


160 


2.54C 


16. 25 


13.75 


.04 


8.0 


0. 9 


25. 6 


161 


2. 54 


16.25 


13.75 


.04 


12.3 


5.8 


57.0 


162 


2.54 


16. 25 


13.75 


.04 


15.8 


4.2 


S3. 5 


163 


2.54 


16.25 


13. 75 


.04 


18.7 


1.8 


88.0 


170 


5.08A 


17.54 


12.46 


.04 


14.0 


7.3 




171 


5.08 


17.54 


12.75 


.04 


16.0 


6. 




172 


5.08 


17. 54 


12. 75 


.04 


19.5 


4.5 




173 


5.08 


17.54 


12.75 


.04 


22.8 


5.2 




174 


5.08 


22.46 


7.54 


.04 


16.2 


9.0 


~"~83.~6 


175 


5.08 


22.46 


7.54 


.04 


19.5 


3.5 


110.0 


176 


5.08 


22.46 


7.54 


.04 


22.4 


3.5 




178 


5.08 


17.54 


12.46 


.04 


13.6 




"""5775 


183 


5.08B 


22.46 


7.54 


.04 


11.9 


~~~4.Y 


41.0 


184 


5.08 


22.46 


7.54 


.04 


10.8 




40.0 


185 


5.08 


22.46 


7.54 


.04 


9.4 


3.3 


25.0 


186 


5.08 


22. 46 


7.54 


.04 


14.2 


6.3 


67.0 


187 


5.08 


22.46 


7.54 


.04 


16.4 


8.3 


89.0 


188 


5.08 


22.46 


7.54 


.04 


19.5 


4.5 


115.0 


189 


5.08 


22.46 


7.54 


.04 


23.3 


2.2 


110.0 


196 


5.08 


17.54 


12.46 


.04 


13.7 




61.0 


197 


5.08 


17.54 


12. 46 


.04 


16.1 




78.0 


198 


5.08 


17.54 


12.46 


.04 


19.6 




90.0 


210 


1.27B 


25 


5 


.04 


8.0 


... .... 


21.6 


211 


1.27 


25 


5 


.04 


5.3 


0.5 


11.3 


212 


1.27 


25 


5 


.04 


9.5 


1.5 


29.2 


213 


1.27 


25 


5 


.04 


10.7 


1.6 


38.2 


214 


1.27 


25 


5 


.04 


12.0 


1.6 


48.8 


215 


1.27B 


25 


5 


.04 


13.8 


1.8 


64. 5 


216 


1.27 


25 


5 


.04 


16.3 


2.0 




217 


1.27 


25 


5 


.04 


L9.6 


1.7 


135. 


218 


1.27 


25 


5 


.04 


22.4 


0.9 


127.0 


219 


1.27 


23. 12 


6.88 


.04 


5. 2 


0.1 


10.0 


220 


1.27 


23.12 


6.98 


.04 


8.0 


0.8 


21.3 


221 


1.27 


23. 12 


6.88 


.04 


9.5 


1.2 


29.5 


222 


1.27 


23.12 


6.88 


.04 


10. 6 


1.3 


36. 


223 


1.27 


23.12 


6.88 


.04 


12.1 


1.5 


47.3 


224 


1.27 


23. 12 


6.88 


.04 


13.8 


1.6 


64.0 


225 


1.27 


23.12 


6.88 


.04 


16. 3 


2.3 


94.0 


226 


1.27 


23.12 


6.88 


.04 


19. 3 


1.5 


146.5 


227 


1. 27 


23. 12 


6.88 


.04 


21.5 


0.9 




228 


L27 


24 


6 


.04 


5.3 


0.2 


10.4 


229 


1.27 


24 


6 


.04 


7.8 


0.9 


19.8 


230 


1.27 


24 


6 


.04 


9.6 


1.3 


29.3 


231 


1.27 


24 


6 


.04 


10. 5 


1.5 


35.1 


232 


1.27 


24 


6 


.04 


11.8 


1.5 


45.0 


233 


1.27 


24 


6 


.04 


13.8 


1.8 


63.0 


234 


1.27 


24 


6 


.04 


16. 1 


2.0 


90.5 


235 


1.27 


24 


6 


.04 


19.2 


1.9 


146. 


236 


1.27 


24 


6 


.04 


22.2 


0.8 


162. 5 


237 


1.27 


20 


10 


.04 


6.3 


0.2 


14.0 


238 


1.27 


20 


10 


.04 


7.9 


0.4 


20.9 


239 


1.27 


20 


10 


.04 


9.5 


0.8 


30.5 


240 


1.27 


20 


10 


.04 


10.7 


1.0 


37.7 


241 


1.27 


20 


10 


.04 


12. 2 


1.2 


49.2 


242 


1.27B 


20 


10 


.04 


13.8 


1.5 


66.0 


243 


1.27 


20 


10 


.04 


16. 2 


2.3 


90.2 


244 


1.27 


20 


10 


.04 


19.5 


1.0 


160.5 


245 

246 


1.27 
1.27 


15. 635 
15. 635 


14. 365 


.04 


6.2 


0.1 




14! 365 


.04 


7.9 


0.2 


~~~22.~6 


247 


1.27 


15. 635 


14. 365 


.04 


9.6 


0.8 


31.6 


248 


1.27 


15.635 


14. 365 


.04 


10.9 


1.3 


41.0 


249 


1.27 


15. 635 


14. 365 


.04 


12.2 


1.9 


53.0 


250 


1.27 


15. 635 


14. 365 


.04 


14.1 


3.3 


77.5 


251 


1.27 


15. 635 


14. 365 


.04 


16.3 


2.0 


126.0 


252 


1.27 


15. 635 


14. 365 


.04 


19.3 


1.0 


138.4 


275 


1. 27A 


20 


10 


.04 


11.8 


1.0 


45.0 


276 


1.27 


15. 635 


14. 365 


.04 


10.8 


1.2 


39.0 


277 


1.27 


15. 635 


14. 365 


.04 


14.1 


3.7 


75.5 



« Cylinders are tabulated by radius (in centimeters) and by shape; A circular 
B tear drop, C lenticular. 
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in which we are interested the equation may be 
written 



where 



a" P N 1 -a 2 pgkN 2 +Apg 2 k 2 =0 

JVi=coth kh coth Jch' + l 
Af 2 =coth M+coth W 

k=2w/\ 

a=kU. 



The equation reduces to 

U 2 N, 






P ' 

where e=Ap/p. There are two solutions of this 
equation corresponding to the plus and minus signs 
for the last term. The minus sign represents the 
case of interest here where waves are produced 
primarily at the interface. Since the quantity 
Nit/N 2 2 is small, we may expand the square root 
obtaining. 

U 2 

TJchL A^ 2 J 



Ap 



gh N 2 l 



For the present experiments Nie/N 2 2 <C0.03, so only 
the first term is taken 



U 2 

=?gh 2irh 

p 



1 



. 2wh , ,, 2wh' 
coth — — )- coth — r— 

A X 



This relation is the general expression for wave 



velocity of a two-layer system. Taking the limit as 

\/h^> oo we obtain 

TJ\ = h' 

Ap , h+h' (a) 

— gh 

p 

for the velocity of long waves. For the condition 
of the experiments, h! = h, we have 

U 2 1 X , , 2irh 

Ap , 4:irh X 

— gh 

p 

which is the basis of the theoretical curves in figures 
6,7, and 8. The agreement is good when U 2 /(Ap/p)gh 
does not exceed 0.3. For larger values of U 2 j 
(Ap/p)gh the wavelength tends to be constant. At 
the present time no explanation can be given for 
this result. 

For the wave heights, neglecting viscous effects, 
we have 

a=a ( U, r, h, s, /3, — g, Lj 

or 

U 2 



:=i 



/ U 2 
\ p 




A plot of a/r against U 2 /(Ap/p)gs is given in figure 9. 
Those points corresponding to cases where the 
cylinder touches the free surface, (s +r)/A' = l, would 
not correlate well with the other points and are 
plotted separately in figure 10. An explanation for 
this will appear in the theoretical development where 



r = 2 54 cm 
> h -- 26.0 cm, CIRCULAR CYLINDER 
O h = 25.0 cm, 
C> h= 24.0 cm, 
O h = 23.0 cm, 
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OX\-- 18.0 cm, 
e h* 25.0 cm, TEAR DROP 
9 h= 20.0 cm, 

<D h= 25.0 cm, LENTICULAR " 
<p h= 20.C cm, 



r = 1.27 cm 

Q h = 20.0 cm, 



B 

B 
B 

r = 

A 
A 



h = 25.0 cm, 
h - 24.0 cm , 
h = 23.1 cm, 
h = 20.0 cm, 
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■ 22.5 cm 
= 22.5 cm 
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Figure 9. Variation of the wave height of internal progressive waves for those cases where the 

cylinder does not touch the surface. 
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Figure 10. Variation of the wave height of internal progressive 
waves for those cases where the cylinder touches the surface. 

successive images of the cylinder in the surface are 
taken. A plot of air against U 2 /(Ap/p)gh is given in 
figure 11. The maximum value of air occurs at 
about U 2 l(Ap/p)gh=0.^. Comparing this with the 
expression for wave velocity U c of long waves it is 
seen that UJU C =0.77, where U m is the velocity 
corresponding to maximum wave height. In this 
respect there appears to be a close agreement with 
Ekman's results of maximum added resistance, as 
large wave height reflects the largeness of the added 
resistance. The data of a/r from runs where the 
cylinder was not in contact with the surface are 
plotted against 



U 2 



Ap 



for 



gs 



JZ! 

Ap 

— qh 
P 



<0.3 



in figure 12. 

From the data on wavelengths and wave heights 
it would be difficult to distinguish any significant 












Figure 11. Variation of the wave height of internal progressive 
waves where the cylinder does not touch the surface. 
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Figure 12. Variation of the wave height of internal progressive 
waves where the cylinder does not touch the surface for 
U*I(Aplp)gh<Q.Z. 

difference between the waves produced by tin 4 
cylinders of various shapes. This may be due to the 
wake being independent of the cylinder shape, al- 
though no observations were made of the wake. 

For sufficiently large values of rJ8 there exists 
a minimum velocity below which progressive waves 
are not produced. Instead, there is a vortex forma- 
tion below the cylinder followed by a smooth un- 
disturbed boundary. It can be seen from figure 12 
that the wave height increases as r/s is increased. 

For large values of rjh the disturbance assumes 
the form of an internal hydraulic jump which will be 
discussed in section 6, 

The effect of the parameter Ljh cannot be evalu- 
ated at this time since L was nearly constant for 
all the runs. 

5. Translation Waves of the Bottom Layer 

For those runs for which h<^<Ch! ', translation 
waves of the form illustrated in figure L3 were formed 
at the interface. In all the runs of this class the 
cylinder touches the free surface. If the cylinder 
were completely below the surface, waves of the 
same form would be produced. The characteristics 
of the intumescence are similar to those of internal 
solitary waves. Letting a be the height of the wave 
crest above the undisturbed level of the interface, 
we have, neglecting viscous effects 

a=a(lJ,r, h,h' , s, (3, — g 7 LJ 



or 




U 



Ap ' h' ti h' P ' h I 

The values of the variables for the runs of this class 
are given in table 2. 

In figure 14 a plot is given of a/r against U 2 / 
(Ap/p)gh. Examining the figure, the maximum value 
of a\r occurs for U 2 /(Aplp)gh=0.52. Comparing this 
with the velocity for long waves, eq (a), the ratio 
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Figure 13. Example of translation waves of the bottom layer produced when h<^<^h'. 
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Figure 14. Variation of the wave height of translation waves of 
the bottom layer. 





Table 2. 


Variables for h< <iY 








—=0.04 for all runs 
p 






Run 


Cylinder" 


h' 


h 


s 


u 


a 






cm 


cm 


cm 


cm/ sec 


cm 


131 


2.54A 


22.5 


10 


20 


9.6 


1.7 


132 


2.54 


22.5 


10 


20 


11.0 


3.8 


133 


2.54 


22.5 


10 


20 


12.7 


7.3 


134 


2.54 


22.5 


10 


20 


14.3 


8.4 


135 


2.54 


22.5 


10 


20 


15.7 


6.0 


139 


2.54 


22.5 


10 


20 


19.8 


2.7 


140 


2.54 


22.5 


10 


20 


22.2 


3.0 


141 


2.54 


22.5 


10 


20 


17.7 


3.5 


152 


2. 54B 


22.5 


10 


20 


14.4 


7.5 


164 


2.54C 


22.5 


10 


20 


14.8 


6.7 


180 


5.08A 


27.62 


7.46 


22.54 


13.7 


14.0 


181 


5.08 


27.62 


7.46 


22.54 


16.2 


8.6 


199 


5.08B 


27.62 


7.46 


22.54 


7.9 


1.8 


200 


5.08 


27.62 


7.46 


22.54 


9.5 


3.0 


201 


5.08 


27.62 


7.46 


22.54 


10.7 


5.0 


207 


5.08 


27.62 


7.46 


22.54 


13.8 


13.0 


208 


5.08 


27.62 


7.46 


22. 54 


15.9 


4.6 


278 


1.27B 


21.27 


10 


20 


10.8 


1.2 


279 


1.27 


21.27 


10 


20 


13.8 


4.0 


281 


1.27 


21.27 


10 


20 


18.6 


2.0 



UJU e has the values 0.87, 0.87, and 0.81 corre- 
sponding to the cylinders r=1.27 cm, 2.54 cm, and 
5.08 cm, respectively. 

For the runs of this class, /*,<</*/, the cylinder 
touches the surface so the value of r/s is small in 
nearly all cases. Thus for low velocities, waves are 
produced and these are of small wave height. For 
these cases progressive waves are produced as in 
section 4. The apparent criterion for transition 
from progressive waves to translation type waves 
is the following: For low velocities progressive waves 
of small amplitude are produced. For larger veloci- 
ties the amplitudes increase until they become sig- 
nificantly large in comparison to the depth h of the 
lower layer for the translation type waves to be 
produced. No effort was made to determine pre- 
cisely these transition velocities. 

6. Hydraulic Jump of the Upper Layer 

For those runs for which A'«A, translation 
waves of the form illustrated in figure 15 were 
formed at the interface. Here again the cylinder 
touches the free surface in all cases. At the start 
of a run progressive waves are formed. As the 
cylinder moves on, the waves grow in amplitude and 
the first crest flattens out because of the proximity 
of the fresh water surface. Letting A be the hori- 
zontal distance between the center of the cylinder 
and the first trough, A increases with the distance 
of travel L. In general, neglecting viscous effects, 



"Cylinders are tabulated by radius (in centimeters) and by shape; A circular 
B tear drop, C lenticular. 



A=A 



(u,r,h',h,^g,s,fi,Ly 



The values of the variables for the runs of this class 
are given in table 3. 

We wish A/L as a function of dimensionless com- 
binations of the variables. The best correlation 
that was obtained is given in figure 16 where A/L is 
plotted against 

■y/r/L. 



Ap 



gh' 



Here again, as in section 5, there are cases, for low 
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Figure 15. Example of internal hydraulic jump of the upper layer produced when h'<<^h. 
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FifjURE 16. Variation of the length of the "plateau" of translation waves of the upper 

layer. 



velocities where progressive waves are produced at 
the interface. There are also intermediate cases 
where the first crest is flattened, but in other respects 
the profile has the characteristics of a progressive 
wave train. There are still further cases when r/h 
is large where the hydraulic jump is produced when 
the two liquids are of equal depth (runs 178, 196, 
and 197 listed in table 3). 

7. Theoretical Development 

In order to see to what degree the observed wave 
heights are in conformity with the laws of hydro- 
dynamics, the case of a cylinder moving with uniform 
horizontal velocity in the upper layer is considered. 
This analysis refers to the oscillatory waves of the 
interface. Corresponding analyses of the translation 
waves of the bottom and upper layers are being 
considered. 



It will be assumed that r is small compared with 
the distance s (see fig. 5). Since the speed of motion 
U is many times smaller than the quantity ->J(j{h-\-h') 
the surface disturbance due to motion of the cylinder 
as it would be in a homogeneous liquid would be un- 
appreciable, and also owing to small Ap/p ratio the 
surface displacement associated with internal waves 
would be negligible. In this case, when the cylinder 
is near the free surface, it will be necessary to take 
the successive images of the cylinder in the surface 
as if it were a plane boundary in order to reduce the 
normal velocity at the surface to zero. For the 
present purpose, since r/s is assumed small, all of 
the image doublets within the cylinder may be 
placed at the center which effectively increases the 
radius of the cylinder by a factor y. The values of 
7 as a function of {h f —s)jr are determined by the 
methods of reference [22] and plotted in figure 17. 
The sharp rise in y for (h'—s)/r=L accounts for the 
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Table 3. Van 


ables for h'<<Ch 






Ap 

— =0.04 for all runs 
p 




Run 


Cylin- 
der a 


h' 


h 


s 


17 


Ai 


u 


A 2 


u 


A 3 


Lz 




119 


2.54A 

2.54 

2.54 

2.54 

2.54 

2.54 

2.54 

2.54 

2.54B 

2.54C 


cm 

7.5 
12.5 
12.5 
12.5 

7.5 

7.5 
7.5 
7.5 
7.5 
7.5 


cm 
25 
20 
20 
20 
25 

25 
25 
25 
25 
25 


cm 
5 

10 
10 
10 
5 

5 
5 
5 
5 

5 


cm/sec 
9.4 
12.6 
13.7 
16.3 
10.8 

12.1 
13.2 
15.5 
13.1 
13.1 


cm 


cm 


cm 


cm 


cm 
73.0 
94.0 
114.0 
115.5 


cm 
867.5 
839.0 
842.0 
819.0 


124 
125 
126 
128 

129 
130 
136 
151 

159 


























91.0 

96.0 
107.0 
136.0 
91.5 
99.0 


517.0 

363.0 

309. 
415.0 
284.5 
310.0 


96.0 

121.5 
151.8 
206. 
131.5 
142.0 


557.5 

572.0 
622.0 
613.0 
553.0 
571.0 


164.0 
187.0 


851.0 
848.0 


183.0 
186.0 


846.0 
839.0 


178 
196 
197 
203 
204 


5.08A 

5.08B 

5.08 

5.08B 

5.08 


17.54 
17.54 
17.54 
12.62 
12.62 


17.54 
17.54 
17.54 
22.46 
22.46 


12.46 
12.46 
12.46 

7.54 
7.54 


13.6 
13.7 
16.1 
14.0 
16.6 


90.0 
82.0 
112.0 
100.0 
131.0 


344.0 
331. 
352.0 
317.5 
332.5 


102.0 
94.0 
142.0 
138.0 


615.0 
582.0 
614. 5 
583.0 


106.0 
111.0 
159.0 
176.0 
208.0 


834.5 
840. 
867.5 
848.0 
854.0 


205 
206 
253 
255 

256 


5.08 

5.08 

1.27B 

1.27 

1.27 


12.62 
12.62 
3.77 
3.77 
3.77 


22.46 
22.46 
27.5 
27.5 

27.5 


7.54 

7.54 

2.5 

2.5 

2.5 


19.5 

22.6 

6.2 

6.7 

7.8 


177.0 


336.0 






328.0 
370.0 
20.5 
17.5 
152.0 


840.0 
824.0 

847.5 
828.5 
841.0 






20.5 
15.5 
78.0 


331.5 
324.0 
322.0 


20.0 
15.5 
115.0 


591.0 
600.5 
551.0 


257 
258 
259 
260 
263 


1.27 
1.27 
1.27 
1.27 
1.27 


3.77 
3.77 
3.77 
3.77 
6.27 


27.5 
27.5 
27.5 
27.5 
25 


2.5 
2.5 
2.5 
2.5 
5 


9.6 
10.6 
11.8 
13.8 

10.5 


68.0 
66.0 
91.0 
131.0 

58.0 


294.5 
292.0 
305.0 
317.5 
303.5 


112.0 
88.0 
126.0 
204.0 
73.5 


573.5 
565.0 
590.0 
609. 
568.5 


152. 
118.0 
156.0 
259.0 
90.0 


837.0 
840.0 
839.0 
840.0 
817.5 




264 
269 
272 
274 


1.27 
1.27 
1.27 
1.27A 


6.27 
8.15 
11.27 
3.77 


25 

23.12 
20 
27.5 


5 

6.88 
10 
2.5 


13.8 
13.9 
13.8 
11.4 


107.0 
92.0 
77.5 

100.2 


302.5 
317.0 
324. 5 
294.3 


148.0 
113.0 
80.5 

147.5 


571.0 
590.0 
575.0 
590.2 


186.0 
131.0 
84.0 
188.5 


865.0 
826.0 
829.5 
857.5 


1 


a Cylinders are tabulated by radius (in centimeters) an 


d by shape; A circular, B tear drop, C lenticular. 








provided x' + xj is negligible there. The origin is 
taken at the undisturbed level of the interface 
directly below the center of the cylinder. The y— 
coordinate is positive upwards and x is positive to 
the right in figure 5. 

If r/(s-\-h) is small the normal component of 
velocity at the bottom from the terms in the bracket 
would vanish. Again X+X! and x'+xj are so 
selected that the normal components of velocity due 
to these terms vanish at the surface and at the 
bottom. Thus assume 
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Figure 17. Variation of radius factor with distance below the 
surface. 


X'= J a! (k) cosh k(y—h') sin kxdk (3) 


increased wave heights when the cylinder touches 
the free surface found in section 4. 

Letting <// be the velocity potential of the upper 
layer and </> that of the lower we write, after super- 
posing a velocity U on the system, 


/> 00 

x[= J a[(k) cosh k{y—h f ) cos kxdk 
and 

r* oo 


<t> f = -Ux 


r\ i m 2 , m 2 i 


+x / +x / 1 


X= a(k) cosh k(y-\-h) sin kxdk (4) 

/.CO 


L ' x 2 +(s-2/)^ tf+yh'-s-yyj 


(1) 


Xi= a x (k) cosh k(g+h) cos kxdk 


and 




<t>=— Ux 


r, | M 2 - (yr) 2 "I 


+x+x x 


where a'(k) and a(k) are to be determined. For the 
equation of the interface, assumed to be steady, we 
put 


L ' x>+{s- y y ' x> + (2h'-s- y yj 


(2) 


r* oo 


Thus the normal velocity at the cylinder is zero 


V= fi(k) cos kxdx, (5) 


2 


50 



and 



t?i=[ji(i) 



Sill 



kxdk. 



The kinematic boundary conditions at the interface 
require that 

^<<+''>=-[£L--[IL- < 6 > 

Equations (1) and (2) may be written in the form 

^ = _Ux-U(yryr[e- k(s - y) +e- Jc{2h, - s - y) ] 

sin kx dk+x'+xl (7) 



and 



<j>=-Ux-U{yr) 2 f 



-k(s-y) 



+ e~ k 



k(2h'-s-y)] 



mikxdk+X+X! (8) 

for s— 2/>0. Differentiating eq (5) with respect to 
x and eq (3), (4), (7), and (8) with respect to y, eq 
(6) yields 



'(k) sinh kh' = —a(k) sinh M 



(9) 



so that 



a'(fc)= 
and 



=--L {^) + (7r) 2 [e-^+e-^»'- s) ]} (10) 



t/ 



-ft(2ft'-s) 



]}• (ii) 



Also 

a[(k) sinh kh f = a x (k) sinh kh=fi 1 (k) 

So the terms xl, x 1? and 7h correspond to a system 

of free waves and will be dropped from the equations. 

If the pressure at the free surface is taken as zero 



and 






provided that p' I p is nearly unity. Letting p'=p 
for y=v and omitting higher order terms in the 
disturbance 

P ' r {gm+U\yr) 2 k[e- ks +e- k{2h '- s) } 
—Uka'(k) cosh kh f ] cos kxdk 

=p r {(jP(k) + U 2 (yr) 2 k[e- ks +e- k(2 ' l '- s) } 

—Uka(k) cosh kh} cos kxdk (14) 



or 



C{U\ir) 2 k[e 



- ks +e- k(2h '- s) \(p'-p) 



+ gP(k)(p'-p)-p'Uka'(k) cosh W 

+ pUka(k) cosh M} cos kxdk=0. (15) 

Introducing «'(&) and a(k) from eqs (10) and T (l 1) 

U 2 (yr) 2 k [«-*•+«-*»'-•>] [(p'-p) 

+ p' coth M/ + p coth feA]+^(i)[sr(p'-p) 

+£/ 2 & P ' coth kh'+U 2 k P coth M]=0 (16) 

or 

«£) = -(7r) 2 

(P'-P) (A-*) 

J 1 )±i L v 

p' l^+ifccoth^l-p [^2-Jfccoth kh] 



[e~ ks +e- 



(17) 




P ' ^p+k coth kh' \—p —v—k coth kh\ 



, [ 6 -^_|_ e -*(2^- S )j CQS kxdk ( 18 ) 



or 

,= -( 7 r)^[- 



2A/-s n 

+ (2ft'- S ) 2 J 



-( 7 r)M 

Jo / 



9 ' \jj2+k coth jfcA'l-p T^~ 2 -fc coth JWil 

X cos kxdk. (19) 

The integral in eq (19) is improper, the denominator 
becoming zero at the value of k for which 

k[ P C3th kh+p' coth kh'] =^j- 2 (p~p'). (20) 

This value, k =2ir/\, corresponds to the wavelength 
of free interfacial waves at the velocity U, as derived 
in section 4. As U increases the wavelength X also 
increases. Taking the limit of eq (20) as k->0 we 
have 



P , P 9 



\P y ( *\ 
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v = — (yry 



} [_x 2 +s 2 



or 



U c = 



l g( P - P ')hh' 
V ph'+p'h 



(21) 



which is the velocity of propagation of long waves. 
For larger values of U there is no real singularity of 
the integrand in eq (19). Taking the principal value 
of the integral in eq (19) for x positive, 



2h'-s 



^0 s -[- g 



-k a {2h'-8)-\ 



sin fax 



-i *r(p'-p)(^2-*o) 

~x 2 +{2h'-s) 2 y P f [coth hh'-hh' csch 2 kji']+p [coth k h- k h csch 2 k G h] 



Jo / 



(p'—p) fL [sin ms+sin (2mh' — ms) ]e~ 



Jo / 



-p)m [cos ms+cos (2m/*/ — ms)]e~ 



M^+m cot m/*/ — p JL ~ m c °t m ^ p'J ?^2+ m cot m ^' ~" P 772 "~ m cot m ^ 



*?=-(y0 : 



{[: 



(22) 



For large values of x only the second term is signif- 
icant. Since rj is an even function of x, eq (5), the 
corresponding expression with x replaced by —x 
holds for negative x. Adding a system of free waves 
to cancel the disturbance far upstream, the result is 



2h'-8 



x 2 +s 2 ' x 2 +(2h 



'-s)*J + p' 



2<jr(p'-p) (Jp-k^j [e- k o s +e-^ m, ~ s) ] sink x 
coth k Q h' — k h' csch 2 k h f }-\- p [coth. k h— k h czch 2 k h] 



+ etc for x>0 



and 



-M> k 



2h'—8 



+s 2 ' x 2 +(2h'-s) 



] 



+ etc for x<^0. 



Thus there is a local disturbance followed by pro- 
gressive waves of length X, where 

ko 

is determined by eq (20). 

8. Comparison With Observation 

Kecalling the conditions of the theory, it is clear 
that it does not apply to the translation waves of 
either section 5 or 6. For a comparison with the 
observed progressive waves, when h=h\ eq (20) 
gives 

(24) 



. h ,t 2wh Ap qh 
4ir - coth — -=— - Z— 
A A p U 



which can be solved numerically to give \/h as a 
function of U 2 /(Ap/p)gh. 

The values are given by the curves in figures 6, 7, 
and 8. The agreement with experiment is good for 



(23) 



U 2 



Ap 



gh 



-<0.3. 



Regarding the wave height a, we have from eq (23) 



h coth koh—kohcsch 2 k h 



(25) 



approximately. Thus a/r is a function of U 2 l(Aplp\)h 
for given values of s/h and r/h. In figure 18 one of the 
theoretical curves is given along with the corre- 
sponding experimental points. For small wave 
heights the observed values are considerably greater 
than the theoretical values. This disparity may be 
due to the wake produced by the cylinder or the 
limitation of the theory to small waves. For the 
smallest waves which could be measured accurately, 
figure 18, the observed values are approximately 
twice as large as predicted by theory. It appears 
that the disparity between observed and theoretical 
wave heights is largely due to the wake produced by 
the cylinder. 
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Figure 18. Variation of wave height of internal oscillatory 
waves; comparison with theory. 



9. Summary 

The disturbances created at the interface of a 
two-layer liquid system by the horizontal motion of 
a cylinder in the upper layer fall into three categories. 
First, when the layers are of equal thickness, in most 
cases a train of progressive oscillatory weaves is pro- 
duced at the interface. Second, when the depth of 
the lower layer is much less than the depth of the 
fresh-water layer, the profile of the interface usually 
resembles that corresponding to a succession of 
positive internal solitary waves. Third, when the 
depth of the lower layer is much greater than the 
depth of the fresh-water layer, in most cases an 
internal hydraulic jump is produced. The excep- 
tional cases correspond to low cylinder velocities. 
The theory for disturbances of small amplitude pre- 
dicts waves of the type observed in the first category, 
but the observed amplitudes are greater by a factor 
of two in those cases where the conditions of the 
theory are most nearly satisfied. This discrepancy 
is ascribed to the effect of the wake behind the 
cylinder. 

For the progressive oscillatory waves there is a 
distinct change in the manner of variation of the 
wave heights corresponding to U 2 /Ap/p)gh = 0.3 
where the wave height reaches a maximum. For 
U 2 /(Ap/p)(jh<^0.3 the wave height a/r depends 
primarily on the quantity U 2 /(Ap/p)gs while for 
U 2 /(Ap/p)(jky>Q.3 the wave height depends on 



U 2 /(Ap/p)gh. Figure 12 shows that the wave 
height increases as r/s is increased. 

The experiments with different shapes of cylinders 
did not yield a significant difference in the disturb- 
ances produced. This may be due to the manner 
in which the wake depends on the cylinder shape. 
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Phase shift effects >n Fabry-Perot interferometry, C. J. 
Koester, ,/. Research NBS 64A, No. 8, May-June 1960. 
A method is demonstrated for utilizing in Fabry-Perot 
interferometry the data on reflection phase shift dispersion 
obtained from fringes of equal chromatic order. Unknown 
wavelengths can be calculated from the Fabry-Perot patterns 
obtained with a large etalon spacing, even without prior 
knowledge of the phase shift of the reflecting surfaces. When 
the theoretical phase shift as a function of wavelength is 
known approximately, then the correct orders of interference 
can be determined for both the Fabry-Perot fringes and 
fringes of equal chromatic order. From the wavelengths of 
the latter the phase shift dispersion can be measured to an 
accuracy of about 10A. The method is especially useful for 
reflectors with large dispersion of phase shift, such as multi- 
layers. Results in the visible spectrum are reported for 
aluminum films and a pair of dielectric 15 layer broad-band 
reflectors. 

Separation and determination of small quantities of aluminum 
in steel, B. B. Bendigo and R. K. Bell, J. Research NBS 64A, 

No. 3, May-June 1960. 

A method is described for determining small amounts of 
aluminum (0.01 to 0.3 percent) in carbon and stainless steels. 
A perchloric-sulfuric acid solution of the steel is eleetrolyzed 
in a mercury cathode cell to remove most of the iron, and a 
chloroform extraction removes elements such as aluminum, 
residual iron, and titanium as cupferrates from a solution 
buffered at pH 3.5. These elements are converted from 
cupferrates to perchloratcs; all except aluminum are then 
extracted as cupferrates with chloroform from 4N hydro- 
chloric acid. Aluminum in the acid solution is determined 
photometrically, using aluminon (ammonium aurintricar- 
boxylate) and a wavelength of approximately 540 mj*. 
Accuracy to within ±0.005 percent aluminum is indicated. 

A symmetric continuous poker model, A. J. Goldman and 
J. J. Stone, J. Research NBS 64B, No. 1, January-March 1960. 
Beginning with Von Neumann, mathematicians concerned 
with the rational analysis of conflict situations have realized 
that investigation of accurate or simplified versions of common 
card games leads to techniques and insights applicable to 
situations of military or economic interest. In the present 
paper, a symmetric poker model one stage more complicated 
than the original Von Neumann game is solved, There are 
two bet levels a, b and an ante of 1 unit (a>6>l) ; no raises 
are permitted. The game has a unique optimal strategy, 
wiiich forbids bluffing on a low hand. The limiting case 
6=1 is shown to yield the Von Neumann model. 

Characteristics of 488 megacycles per second radio signal 
reflected from the moon, B. C. Blevis and J. H. Chapman, 

J. Research NBS 64D, No. 4, July-August 1960. 
Radio signals at 488 Mc, received after reflection from the 
moon, have been studied using a continuous wave bistatic 
radar system located near Ottawa, Canada. These experi- 
ments were carried out during 1957 and consisted of monitor- 
ing the signals received on two orthogonal dipoles mounted 
at the focus of a 28-foot parabolic tracking antenna. The 
total mean received signal yielded an effective radar cross 
section of the moon at this frequency of 0.05 of the projected 
area. The libration fading as observed in the two orthogonal 
receiver channels was in sychronism, with a cross-correlation 
coefficient of 0.89. It was established that the signal is not 
depolarized significantly on reflection at the surface of the 
moon, or in passing through the ionosphere. Rotation of the 
plane of polarization of the radio wave in the double passage 
through the ionosphere was observed to change by nearly 



180 degrees over a 6-hour period during quiet ionospheric 
conditions. 

The use of polarization fading of satellite signals to study the 
electron content and irregularities in the ionosphere, C. G. 

Little and R. S. Lawrence, J. Research NBS 64D, No. .' tr 
July-August 1960. 

A procedure is described for using the Faraday-rotation 
fading of a satellite radio signal to measure the ionospheric 
electron content per unit column up to the height of the satellite. 
Full account is taken of ionospheric refraction by using the 
collision-free form of the Appleton-Hartree equation. Re- 
sults based on observations of 1958 Delta 2 are presented. 
The sub-satellite electron contents derived throughout the 
satellite passes for heights below the F-peak compare well 
with values deduced from simultaneous ionograms. Iono- 
spheric irregularities, having lateral dimensions of a few 
hundred kilometers, and fractional deviations in sub-satellite 
electron content of about 0.02 have been detected. Our 
observations suggest that satellite polarization studies ofTer 
important advantages over other methods of investigating 
these irregularities. 

Determination of piezoelectric properties as a function of 
pressure and temperature, J. E. McKinney and C. S. Bowyer, 
J. Acoust. Soc. Am. 32, No. 1, 56 {I960). 
Piezoelectric data on a mixed titanate system (82% BaTi() 3 , 
9.1% CaTiO.3, 3.6% PbTi() 3 , and 4.4% Ti0 2 ) was obtained 
from an apparatus intended to measure the dynamic com- 
pressibility of materials. The calibration constant of the 
apparatus involves the piezoelectric and dielectric constants 
of the ceramics used as transducers. The method is dis- 
cussed and an operating equation for the apparatus derived. 
Piezoelectric and dielectric constants have been measured 
over a fairly extensive temperature and pressure range and an 
apparent phase transition observed near 5° C. This transi- 
tion corresponds to that found in pure barium titanate, and 
might be due to imperfect mixing of the titanates that Leaves 
a small amount of pure barium titanate present. The method 
described appears more sensitive for the detection of change's 
in crystal symmetry in such materials than either X-ray or 
dielectric measurements. 

Determination of the recording performance of a tape from its 
magnetic properties, E. I). Daniel and I. Levine, J. Aeons/. 
Soc. Am. 32, No. 2, 258 {I960). 

This is the second of two papers. The first (I) dealt with the 
general magnetic properties of recording tape. The present 
paper gives an analysis of hf biased recording based upon the 
anhysteretic properties of the tape as specified in I. Calcu- 
lations of performance are compared with the results of abso- 
lute measurements carried out on a conventional recording- 
machine, and on large scale models. Satisfactory agreement 
is obtained when the theory is corrected to take into account a 
loss attributed largely to a dynamic form of demagnetization 
which takes place during recording. The corrected theory 
can be used to compare tape performance under a variety of 
operating conditions without reference to a particular test 
recorder. Attention is confined to recorded wavelengths 
which are large compared with coating thickness. 

Gages for measuring the thickness of chromium on the in- 
ternal surface of small-bore tubes, V. A. Lamb and P. A. 
Krasley, Plating 47, No. 2, 176 {I960). 

Two gages of the eddy-current type for measuring the thick- 
ness of chromium coatings on the internal surface of caliber 
.30 and caliber .50 gun barrels have been developed. Meas- 
urements are nondestructive and are localized to a small area, 
less than x /% inch diameter. The gages may readily be modi- 
fied to be generally applicable to the measurement of the 
thickness of various coatings on the internal surfaces of metal- 
lic tubes of the same or larger diameters. 
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Experimental and theoretical investigation of the magnetic 
properties of iron oxide recording tape, E. D. Daniel and I. 
Levine, /. Acoust. Soc. Am. 32, No. 1, 1 (1060). 

Part A gives the results of remanent magnetization tests made 
under ordinary and anhysteretic conditions, and shows that 
the major anhysteretic properties of a recording tape can be 
expressed in terms of three easily measured constants. The 
design of the test equipment is discussed and test results are 
listed for thirteen representative types of tape. Part B re- 
views some of the theories of fine particle magnets that can be 
applied to recording tape, and gives an extensive treatment 
of remanent magnetization based upon the Preisach diagram. 
Some aspects of the Preisach diagram treatment may be of 
interest to workers outside the magnetic recording field. The 
anhysteretic properties are important in hf biased recording 
and a second paper describes how these properties can be used 
to predict certain recording performance characteristics. 

Efficient harmonic generation, G. F. Montgomery, Proc. IRE 
48, No. 2, 251 (1960). 

An harmonic generator using a rectifier multiplier followed by 
an amplifier is more efficient than its equivalent class-C multi- 
plier. A proposed rectifier-transistor harmonic generator is 
capable of high efficiency that is nearly independent of har- 
monic number. 

Precise measurement of heat of combustion with a bomb 
calorimeter, R. S. Jessup, NBS Monograph 7 (1960). 

This Monograph gives detailed descriptions of apparatus and 
methods which are used at the National Bureau of Standards 
for precise determinations of heats of combustion of liquid 
hydrocarbon fuels. Numerical examples are given of methods 
of calculating results of measurements from observed data. 
The technique of making and filling glass bulbs to contain 
samples of volatile liquid fuels is described. 
The accuracy of the methods described is about 0.1 percent. 
This is intermediate between the accuracy of 0.01 br 0.02 
percent attained in certain measurements on pure com- 
pounds, and the accuracy of several tenths of one percent 
obtainable with published standard procedures for measure- 
ments on fuels. 

Conductive flooring for hospital operating rooms, T. H. Boone, 
F. L. Hermach, E. H. MacArthur, and R. C. McAuliff, 
NBS Monograph 11 (1960). 

Characteristics and performance of available types of con- 
ductive flooring materials were investigated in the laboratory. 
The study showed that the electrodes and instruments 
used to measure the floor greatly affected the measured 
resistance, but that the method specified by the National 
Fire Protection Association for measuring the electrical 
resistance reasonably simulated the conditions under which 
a floor functions in reducing electrostatic hazards. The 
physical, chemical, and serviceability characteristics of 
conductive floorings investigated showed results comparable 
with those of nonconductive flooring of the same type. 
Consequently, with some limitations, an architect may base 
his choice of a conductive flooring material on his knowledge 
of the behavior of similar nonconductive materials. 

Tabulation of data on receiving tubes, C. P. Marsden, W. J. 

Keery, and J. K. Moffitt, NBS Handb. 68 (1959). 
A tabulation of Receiving Type Electron Tubes with some 
characteristics of each type has been prepared in the form 
of two major listings, a Numerical Listing in which the tubes 
are arranged by type number, and a Characteristic Listing 
in which the tubes are arranged by tube type and further 
ordered on the basis of one or two important parameters. 
The tabulation is accompanied by a listing of similar tube 
types and basing diagrams for the listed tubes. 

Are life testing procedures robust? M. Zelen and M. C. 
Dannemiller, Proc. 6th Natl. Symp. Reliability and Quality 
Control in Electronics, Jan. 11-18, I960, Inst. Radio Engrs. 
Inc., p. 185 (1960). 

Almost all the statistical procedures in current use for evaluat- 
ing the reliability of components or equipments rest on the 



assumption that the failure times follow the exponential 
distribution. In practical situations one rarely has enough 
data to determine whether failure times are actually exponen- 
tial. This paper studies the behavior of several statistical 
life testing procedures based on the exponential failure law 
if the true failure law is the Weibull distribution. It is found 
that these statistical techniques, which are widely used, are 
very sensitive to departures from initial assumptions. Apply- 
ing these techniques to life test data when the exponential 
failure law is not satisfied may result in substantially increas- 
ing the probability of accepting components or equipments 
having poor reliability. 

Improved NBS abrasive jet method for measuring abrasion 
resistance of coatings, A. G. Roberts, ASTM Bull. No. 244, 
48 (TP52) (1960). 

A description is given of the changes made in the NBS 
Abrasive Jet Method since publication of the original article 
describing this rapid, reproducible, and versatile means for 
measuring the abrasion resistance of organic and other 
coatings. Improvements include use of a free-flowing abra- 
sive, a redesigned abrasive sifting chamber to minimize 
turbulence and interference effects, better control of abrasive 
flow, addition of an electronic contact device for setting 
nozzle-to-coating distance, and automatic timing. Precision 
of measurement has been significantly increased by use of a 
reference panel. 

Flame-spread measurements by the radiant panel flame- 
spread method, D. Gross, Forest Products J. X, No. 1, 33 

(1960). 

The results of a limited number of tests for evaluation of the 
flame spread properties of untreated and of fire-retardant 
coated wood-base materials are presented. It is shown by 
use of the radiant panel flame spread test method that 
variations in moisture content and mean density can have 
a significant effect upon the flame spread index of untreated 
wood. For the same mean density and moisture content, it 
was observed that variations in the surface grain structure 
of untreated wood had an appreciable effect upon the flame 
spread index. The flame spread indices of the fire-retardant 
coated assemblies were considerably lower than those obtained 
with uncoated or conventionally coated materials although 
significant differences between coated assemblies were 
observed. A test method such as this can be used to ad- 
vantage for the development and evaluation of effective and 
practical fire-retardant coatings for wood-base materials. 

Report of the 44th National Conference on Weights and 
Measures 1959, NBS Misc. Publ. 228 (1959). 

This is a report of the proceedings (edited) of the Forty- 
fourth National Conference on Weights and Measures, spon- 
sored by NBS, held in Washington, D.C., June 8-12, 1959, 
and attended by State, county, and city weights and meas- 
ures officials, and representatives of the Federal Government, 
business, industry, railroads, and associations. 

What price accurate test methods? A. T. McPherson, 
ASTM Bull. (ACR Notes Column) No. 244, 7 (1960). 
Inaccurate test methods used in specifications frequently 
lead to difficulty, misunderstandings, and unnecessary rejec- 
tions because much of the tolerance that is allowed in dimen- 
sions, composition, or other characteristics of a product may 
be taken up by the uncertainty in testing. Accurate methods 
can pay dividends by lowering the rejection rate, avoid- 
ing referee tests and other sources of delay, and affording 
greater assurance that products which are found by test to 
comply with the specification will actually serve the intended 
purpose. A further saving may be effected with many prod- 
ucts by widening the tolerances and thereby reducing the 
cost of production. The development of accurate test meth- 
ods, however, may require long and expensive research. 
Before undertaking to develop any specific new test method, 
a critical survey should be made to determine whether the 
anticipated savings and other benefits will exceed the prob- 
able cost. Such a survey will involve both economic and 
technical factors and may well make use of methods and 
techniques which are employed in market research. 
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Many changes reflected in new dry cell standard, W. J. 

Hamer, Mag. of Standards 31, No. 3, 81 {I960). 
The seventh edition of the American Standard Specification 
for Dry Cells and Batteries was published in December 1959 
as Handbook 71 of the National Bureau of Standards. It 
supersedes the previous specification which \v:is published in 
1954 as Circular 559 of the National Bureau of Standards. 
The new specification includes many changes. Cell dimen- 
sions are given in both the English and metric systems: the 
metric system was added in line with the trend toward its 
increased use in commerce and industry. Cell and bat- 
tery designations recommended by TC-35 Committee on 
Primary Cells and Batteries of the International Electro- 
technical Commission was incorporated in keeping with 
attempts to secure international cooperation in the standard- 
ization of primary batteries. The most notable change in 
the new specification is the inclusion, for the first time, of 
cells and batteries for use with transistorized circuits. Other 
minor changes include a new test for lantern batteries used 
in alarm circuits, some new terminal designs, and the drop- 
ping of specifications for batteries to be used with " carbon- 
type" hearing aids since these types of instruments are 
rapidly becoming obsolete. 

Water penetration testing machine for sole leather, T. J. 

Carter, /. Am. Leather Chem. Assoc. LV, No. 8, 139 {I960). 
A test to measure the water resistance of sole leather has been 
developed. A leather specimen is flexed at a frequency of 
30 cycles per minute with one side exposed to water until 
penetration occurs. The specimen is clamped so that the 
angle of bend is approximately 90°. The results show that 
flexing increases the rate of water absorption as well as the 
rate of water penetration compared to the rates under static 
conditions. Studies were made on the influence of thick- 
ness of specimen, location on the sole leather bend, and com- 
pression of the specimen on dynamic water penetration and 
water absorption. The logarithm of the number of flexes 
required for penetration is a direct function of the thickness 
of the leather. The rates of water penetration and water 
absorption depend on the location of the specimen on the 
sole leather bend. Compression of the impregnated leather 
decreases the rate of water penetration and the rate of water 
absorption. 

Bending and stretching of corrugated diaphragms, R. F. 

Dressier, Trans. ASM E SID, No. 4, 651 {1959). 
Solutions of the exact linear elastic shell equations for all 
stresses and displacements in a typical corrugated diaphragm 
are presented for three cases over a 9 to 1 range of thicknesses. 
Results were obtained by numerical integration in an elec- 
tronic digital computer. The effect of thickness variation is 
discussed with respect to both stresses and resultants, and 
peak values needed for design purposes are presented. Cir- 
cumferential and meridional stresses are found to be equally 
important throughout the thickness range analyzed. Bending 
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Triaxial tension at the head of a rapidly running crack in a 
plate, J. M. Frankland, Trans. ASME, Ser. E 26, No. 4, 570 

(1959). 

Two-dimensional solutions have been found for the state of 
stress in a plate loaded in tension at right angles to a crack. 
One of these concerns the stationary crack and the other an 
advancing crack. No three-dimensional solutions are avail- 
able. It is shown, however, that it is highly probable that 
cracks propagating at high speeds develop an appreciable 
triaxial tension near the head of the crack. This is believed 
to be significant with regard to the observed brittle behavior 
of plates under such circumstances. 

Transistor P-A amplifier, G. F. Montgomery and F. R. 
Bretemps, Electronic Ind. 19, No. 1, 196 {I960). 
A portable public-address amplifier using transistors and 
powered by flashlight cells develops an electrical power of 
about 0.75 watt and is useful for moderately large audiences. 
The complete unit consists of the amplifier, two loudspeakers, 
batteries, and a dynamic microphone and is contained in an 
attache case. Battery life is about 20 hours. 
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